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Abstract 

We interpret the discrepancy between semiempirical nuclear mass for- 

^T^ ■ mulas and actual nuclear masses in terms of the residual interaction. We 

r-{ ' show that correlations exist among all binding energies and all separation 

rH , energies throughout the valley of stability. We relate our approach to 

chaotic motion in nuclei. 
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C^ ' Motivation. The masses of atomic nuclei are keys to the understanding of 

many physical and astrophysical processes. For this reason it is important to 
construct reliable theoretical models for the values of the nuclear masses or, 
equivalently, of the binding energy B{A) as function of mass number A. (For 
simplicity of notation we suppress the additional dependence on proton number 
Z). This function is also needed to predict the masses of (as yet) unknown 
nuclei. 

The standard approach to a global modelling of the function B{A) starts 
out from the liquid-drop model for the nucleus and considers in addition shell 
corrections as well as corrections due to pairing. The resulting semi-empirical 
mass formula contains about 30 parameters and is fitted to a large number 
of data. Years of painstaking work have culminated in a best fit P^ which 
reproduces the data points very well albeit not exactly. The overall difference 
between that best fit and the actual data is of the order of 0.5 Mev, with 



a tendency to become smaller for heavier nuclei. That figure, although very 
small in comparison with the total binding energy, is not negligible. Other 
approaches [21 E] lead to similar differences. 

Bohigas and Leboeuf ^ have suggested that the discrepancy can be at- 
tributed to the fact that the nuclear dynamics is partly chaotic. Using the 
semiclassical approximation and Berry's result |S] for the form factor of the 
two-point correlation function, these authors have estimated the contribution 
to the variance of the binding energy due to chaotic single-particle motion. For 
the square root of the variance, the estimate yields 

Like the shell-correction employed in the best-fit procedure of Ref. PP , Eq. ^ 
is derived in the framework of the single-particle picture and is, thus, based 
upon a mean-field description. However, the authors of Ref. fj state that this 
equation "presumably gives an estimate of the mass fiuctuations from neglected 
many-body effects" . When plotted versus mass number, the result (^ is in 
good agreement with the observed discrepancy except for the lightest nuclei, 
see Fig. 1 of Ref. (H. 

The work of Ref. ^ has led to a lively discussion in the community. In 
this letter, we wish to address some of the points that have come up in the 
discussion. Our point of view differs from that of Ref. ^ in that we interpret 
the discrepancy in terms of many-body effects due to the residual interaction of 
the shell model and explore the consequences of such an interpretation. We are 
primarily interested in the conceptual implications of such an approach. Many- 
body effects have been addressed in a number of publications, see, for instance, 
Refs. 000 El- We discuss some of these works below. 

The Two-Body Random Ensemble. At the outset, we have to specify the 
residual interaction. We assume that the residual interaction is of two-body 
form. We do so only to have a specific model within which we can work. We 
are aware of the fact that there is evidence for weak three-body interactions in 
nuclei jJOl ■ It will be seen that allowing for three-body forces will not affect our 
conclusions. 

Within every major shell, the residual two-body interaction is specified in 
terms of a set of two-body matrix elements. All possible two-body interac- 
tions are taken into account by considering the two-body matrix elements as 
Gaussian-distributed random variables. This defines the two-body random en- 
semble (TBRE) of the shell model |nilI21- The TBRE has the advantage of 
leading to generic statements (valid for almost all two-body interactions). This 
is useful since not all of the two-body matrix elements of the residual interac- 
tion of the shell model are precisely known either empirically or theoretically, 
especially in medium- weight and heavy nuclei, and for higher-lying shells in 
any nucleus. Moreover, the use of the TBRE allows us to estimate variances 
of observables and, thus, the range within which values of the observables will 
vary as we vary the parameters of the residual interaction or, equivalently, to 



estimate the uncertainty of the observables due to incomplete knowledge of the 
residual interaction. 

If one applies the TBRE to the calculation of binding energies (more pre- 
cisely: to estimate corrections to the semi-empirical mass formula as employed 
in Ref. ^), the binding energies will vary from one realization of the TBRE to 
the next. Technically speaking, the binding energies become random variables. 
Again, their variances express our uncertainty about the residual interaction. 
We use these variances as a measure of the deviations from empirical values of 
the binding energies that are expected to occur in semi-empirical mass formulas 
which neglect parts or all of the residual interaction. (The actual formula used 
in Ref. ,Ij does not disregard the residual interaction altogether but does take 
into account the pairing force and the symmetry energy. The version of the 
TBRE used here must take this fact into account by keeping fixed the corre- 
sponding two-body matrix elements, and by considering as random only the 
remainder). 

Binding Energies and Separation Energies. We discuss the application of 
the TBRE to binding energies. We first recall the connection between nuclear 
binding energies and nucleon separation energies. We denote the separation 
energy of the least-bound nucleon in nucleus A (which is either a neutron or 
a proton) by <S'(^), again omitting for simplicity the dependence of S on Z. 
We can completely decompose the nucleus A into its constituent nucleons by 
removing one nucleon after the other. Thus we have the identity 

A 

B{A)=Y,SU). (2) 

Likewise, since after removal of k nucleons we reach the nucleus with mass 
number A — k (with k integer), we also have 

A 

B{A)^B{A-k)+ Y. ^O')- (3) 

j=A-k+l 

Given the initial nucleus with mass number A and proton number Z , the S{j) 
in Eqs. Q and Q are uniquely specified by the requirement that we always 
remove the least-bound nucleon. With this understanding, Eqs. (0) and Q hold 
for all nuclei (and not only for those in the bottom of the valley of stability) 
and for all values of k. More equations of the same type are obtained when we 
allow for the removal of the least-bound proton or least-bound neutron when 
that is not the least-bound nucleon. Hence, equations of the type of Eqs. ^ 
and (0 constitute a dense network of identities that hold throughout the mass 
valley. 

It was stated above that as we apply the TBRE to calculate corrections to 
the binding energies predicted by the model of Ref. J4, such binding energies 
become random variables. The same holds for the nucleon separation energies. 
For clarity these are denoted by the letters B and S, respectively. The random 



variables B{A) and S{j) must obey equations of the type ||2Jl and Q, since 
by the very definition of the binding energy and the separation energy, every 
reahzation of B{A) and S{j) in the framework of the TBRE must obey these 
equations. Thus, 

A 

B{A) = J2s{j), (4) 

A 

B{A)= J2 S{j) + B{A-k). (5) 

j=A-k+l 

Correlations. We recall that the deviations of the data from the empirical 

mass formula decrease with increasing values of mass number A, see Eq. ^. 

We now show that this fact combined with Eqs. Q) and |(SJ) implies strong 

correlations between the random variables B{A) and S{A) for all values of A. 

We recall that stochastic contributions to these variables arise from the TBRE. 

If we assume that all variables B{A) and S{A) are uncorrelated, it follows from 

Eq. gl) that 

A 

Var(e(A))=^Var(5(j)). (6) 

The variance is, by definition, positive semidefinite. Therefore, Eq. (|SJ| implies 
that Var(S(^)) is a non-decreasing function of A. (That same conclusion was 
drawn already in Ref. ^] albeit in the framework of a specific model. Here 
we establish it in full generality). This result is in striking contrast to the 
monotonic decrease of the discrepancy with increasing mass number displayed 
by the data and by Eq. (^). We conclude that there exist strong correlations 
between binding energies and separation energies. 

Allowing for such correlations to exist, we arrive at a modified form of 
Eq. (0). With the covariance of two random variables X and Y defined as 
Cov(xr) = Xy - X y, we have 

A A 

Var(B(A))=^Var(5(j))+ ^ Cov(5(z)5(j)) . (7) 

In order to attain the observed monotonic decrease of the discrepancy with 
mass number, the sum of the covariances must be negative semidefinite for all 
values of A and cannot vanish identically. This definitely implies the existence 
of correlations among the <S(j)'s. Given this fact, we expect correlations to exist 
also among binding energies pertaining to different values of A. Indeed, Eq. Q) 
implies 

A-k A A-k 

Cov{B{A)B{A - fc)) = ^ Var(5(j)) + E E (^ " Kn)Coy{S{n)S{j2)) • 

(8) 



This expression is, in general, different from zero. In contrast to the case of 
the separation energies we have, however, no firm empirical evidence in the 
present framework that these correlations do indeed exist. They may vanish 
accidentally for all values of A. Even in this case we expect correlations of 
higher order than the second to differ from zero. Indeed, it takes a very special 
form of the probability density of random variables connected by linear relations 
like Eq. @ to be uncorrelated. We conclude that within the framework of the 
TBRE, all binding energies and all separation energies throughout the mass 
valley are correlated. 

Discussion. From the point of view of nuclear physics, this result is not 
terribly surprising. It mirrors the fact that a change of the residual interaction 
within a major shell will affect the binding energies and the separation energies 
of all nuclei pertaining to that shell. As we cross the boundary in mass number 
separating major shells, the dominant two-body matrix elements of the residual 
interaction change, and we might expect a weakening of the correlations. Our 
result shows no trace of such an effect, and this is probably its most surprising 
feature from the nuclear physics point of view. 

Unfortunately, we cannot draw any conclusions about the range of the cor- 
relations. Still, the following model is instructive. We assume that the cor- 
relations between the separation energies are minimal and restricted to near- 
est neighbors only so that Cov{S{ji)S{J2)) oc Sj^j2±i- For the range of the 
correlations, this is a worst-case scenario. Our assumption is consistent with 
Eq. (UJ provided all these covariances are negative definite. By a suitable choice 
of the values of the covariances, it is possible to account for the monotonic 
decrease of a with A as postulated in Eq. (^. However, even that model im- 
plies long-range correlations between the binding energies. Indeed, under the 
assumptions of the model, the right-hand side of Eq. IjSJl becomes equal to 
Var(B(A - k)) + Cov{S{A - fc + 1)S{A - k). Thus, the covariances of B{A - k) 
with B{A) have for all values of A the same values! We have not been able to 
design a model which would satisfy Eqs. iQ and © with short-range corre- 
lations between all 5's and all B's. From the point of view of the TBRE, we 
expect the correlations to extend over a range of mass numbers given by the 
range of major shells, see Ref. pH. 

Statistical fluctuations of nuclear spectra are commonly described not in the 
framework of the TBRE but rather in the framework of the Gaussian orthogo- 
nal ensemble of random matrices (GOE). From the point of view of canonical 
random-matrix theory, the existence of long-range correlations is unexpected. 
Indeed, at present there exists no conceptual framework within the GOE in 
which the existence of such correlations could be accommodated. Thus, the 
data point to a further limitation of canonical random-matrix theory. Earlier 
work on sd-shell nuclei |l3j has revealed such limitations in a more restricted 
framework by proving the existence of correlations between levels in the same 
nucleus carrying different quantum numbers and between levels in different nu- 
clei belonging to the same major shell. 

What is the relation between our point of view - which blames the discrep- 
ancy between data and the nuclear mass formula on the residual interaction - 



and that taken in Ref. ^ - which puts the blame on order-chaos coexistence? 
We identify chaos in nuclei with the existence of spectral fluctuations of the 
Wigner-Dyson type. This statement applies to levels of fixed spin and parity 
within the same nucleus. In the following three paragraphs we present argu- 
ments to the effect that such spectral fluctuations are always due to the residual 
interaction. Given this fact, we conclude that the two points of view are not 
fundamentally different especially in view of the remark ^ that Eq. |^ "pre- 
sumably gives an estimate of the mass fluctuations from neglected many-body 
effects". 

To justify our view that the residual interaction is the agent responsible for 
chaos in nuclei, we consider first spherical nuclei. Here the single-particle mo- 
tion is regular and does not lead to level sequences which obey Dyson-Mehta 
statistics. This changes when the residual interaction is included. Extensive nu- 
merical studies by many authors, mainly of nuclei in the sd-shell, have shown 
that in this case, nearest-neighbor spacing distribution and Mehta-Dyson statis- 
tic do agree with GOE predictions whenever the number of nucleons in the shell 
exceeds three or four. This is true provided the matrix elements of the residual 
interaction are sufficiently strong to mix configurations pertaining to different 
single-particle energies within the same major shell. For a review, see Ref. |15) . 
Recent theoretical studies of the TBRE confirm this picture for the sd-shell ^Hl 
and show that almost all two-body residual interactions cause chaos. General- 
izing this argument we are led to expect that in every major spherical shell, a 
sufficiently strong residual interaction generically produces chaos. (This state- 
ment is not at variance with the existence of regular spectral features near the 
ground state or like the giant resonances. It relates to spectral statistics, i.e., to 
the joint properties of a large number of levels). We add in parenthesis that the 
Skyrme-Hartree-Fock approach ^7] or other similar approaches take account 
of many-body effects in terms of an effective single-particle density functional. 
In this approach the line separating single-particle and many-particle effects 
becomes blurred. 

The situation is somewhat different in deformed nuclei. For sufficiently 
strong deformations, the single-particle motion alone may be chaotic, and the 
single-particle spectrum may obey Wigner-Dyson statistics. But strong de- 
formations occur in the middle between large major shells, so we deal with 
many valence nucleons. And when we put several or many nucleons into these 
deformed single-particle states and observe the exclusion principle, we gener- 
ate a spectrum which by construction is a superposition of several or many 
GOE spectra. It is known that such superpositions tend rapidly (with increas- 
ing nucleon number) towards the Poisson spectrum which is characteristic of 
integrable motion. This suggests that also in deformed nuclei, the residual in- 
teraction is needed to mix the configurations and to generate spectral properties 
of the Wigner-Dyson type. We are not aware, however, of studies to show under 
which conditions sufficiently strong mixing is attained. 

A different view of chaos in nuclei is offered by the Interacting Boson Model 
which covers the case of deformed nuclei but also that of nuclei with vibrational 
spectra etc. That is an effective model which describes the low-energy spec- 



tra of nuclei with many nucleons outside of closed shells, i.e., of nuclei where 
shell-model calculations are ineffective, in terms of several interacting bosons of 
angular momenta zero and two. Without any residual interaction between nu- 
cleons, the model would not exist. Studies of chaos within this model JHl have 
shown that chaos is generic except near those parameter values where dynamical 
symmetries persist. 

The available evidence then suggests that except for very special cases where 
symmetries dominate, it is legitimate to say that almost all forms of the resid- 
ual interaction cause chaos, and that chaos would not exist without the resid- 
ual interaction. This view is also supported by the empirical evidence, i.e., 
the agreement of the spectral fluctuations of the nuclear data ensemble with 
Wigner-Dyson statistics ^^l j and the available evidence at lower excitation en- 
ergies, see Ref. j^O] and references therein. 

This view sheds new light on attempts to reduce the discrepancy between 
mass formulas and data by taking into account the residual interaction in some 
way, either by explicit diagonalization, or by using relations between binding en- 
ergies of the Garvey-Kelson type or some other local (in A) information. Several 
successful such attempts have been published, see, for instance, Refs. [H1[3|S1|U]. 
Inasmuch as the residual interaction produces chaos, this success does not con- 
tradict the interpretation of the discrepancy in terms of chaotic nuclear motion. 
The same remark would apply to versions of the energy-density functional the- 
ory where the residual interaction is incorporated in the single-particle Hamil- 
tonian. 

In some papers |21l |H] chaos in nuclei is interpreted as limiting the accuracy 
with which spectra or binding energies in nuclei can be calculated. It is claimed 
that the error estimated in Eq. (^ constitutes a bound on the possible accuracy 
of dynamical calculations of nuclear masses, and any reduction of this bound 
is then taken as proof that the discrepancy does not have its origin in chaos. 
We disagree with this point of view. The spectra of nuclei in the sd-shell, for 
instance, do display chaos but can be calculated without any limitation on the 
numerical accuracy. This point is discussed in Ref. |16| . More generally, it 
was shown in Ref. [22j that while classical chaos does indeed limit dynamical 
calculations to a time scale given by the Lyapunov coefficient, this is not true 
in general for quantum calculations. 

In a recent paper 1281 it was shown that correlations among binding energies 
can also be accounted for in the framework of the mean-field approach. So far, 
this work is restricted to neighboring isotopes differing by a few units in mass 
number. 

In summary, we have interpreted the discrepancy between the nuclear mass 
formula and nuclear masses in terms of effects due to the residual interaction. 
Using the two-body random ensemble, we have shown that correlations exist 
between all binding energies and all separation energies throughout the valley of 
stability. To the best of our knowledge, this is the first evidence for such wide- 
ranging correlations in nuclei. We have argued that chaos in nuclei is generically 
caused by the residual interaction and does not exist without it. Therefore, we 
believe that there is no substantial difference between our view and the one taken 



in Ref. ^. There are successful attempts to reduce the discrepancy by taking 
into account parts of the residual interaction not contained in the semiempirical 
mass formula in theoretical calculations. Such attempts do not contradict the 
interpretation of the discrepancy displayed in Ref. 4 in terms of chaotic nuclear 
motion. 

The authors are grateful to T. Papenbrock for helpful discussions. 
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